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Abstract. We study the threshold > 2 contact process on a homogeneous tree Tf, of 
degree k = b + 1, with infection parameter A > and started from a product measure 
with density p. The corresponding mean-field model displays a discontinuous transition 
at a critical point \^^{k,9) and for A > X^^{k,9) it survives iff p > {k,6, \), where 
this critical density satisfies < {k,6, X) < 1, Yimx^^ p^^ {k,9, X) = 0. For large 
b, we show that the process on T;, has a qualitatively similar behavior when A is small, 
including the behavior at and close to the critical point Ac(Tb,0). In contrast, for large 
A the behavior of the process on Ti, is qualitatively distinct from that of the mean-field 
model in that the critical density has pc(Tb,0,oo) := lim;\^oo pc(Tf,, 0, A) > 0. We also 
show that limf,^oo ^c(T(,, 0) = where 1 < $2 < $3 < limg^oo = 00, and 
0<liminffe^^ 6^/(^-1) Pc(Tfe,6»,cx)) < limsupf,^^ pc(Tb, 0, cxd) < cx). 



1. Introduction and results. 

1.1 Preliminaries. Threshold contact processes form a natural class of interacting 
particle systems (see, e.g., [Ligl] for background to the area). They are most naturally 
defined on a graph or oriented graph G = {V,E). In both cases V is an arbitrary 
countable set, whose elements are called vertices or sites. When G is a graph, E C 
{{v,u} : v,u E V} is its set of edges, also called bonds. When G is an oriented graph, 
E' C y X y is its set of oriented edges, also called oriented bonds. We denote the influence 
neighborhood of f G ^ in G by 
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{u E V : {v, u} G E}, if G is a graph, 

{u E V : {v,u) E E}, if G is an oriented graph. 
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The degree of the site v eV is the cardinahty of Mg,v The threshold 6 contact process 
on G is now defined as the Markov process on {0, 1}^ with flip rates at v e F at time 
t > given by 

• 1 flips to at rate 1. 

• flips to 1 at rate A in case there are at least 9 sites of in state 1 at time 
t, and at rate otherwise. 

The parameter A > is called the infection rate. The state of the process at each site at 
each time is called the spin at that site at that time. A spin is interpreted as a vacant 
or healthy site, while a spin 1 is interpreted as an occupied or infected site. It is well 
known (see, e.g.. Chapter I of [Ligl]) that such rates define a unique Markov process. 
Note that the fiip rates above are attractive (see Chapter III of [Ligl]), a property that 
has many consequences. 

When ^ = 1, the threshold contact process is easier to analyze, among other reasons 
because it has an additive dual process. (In this dual process infected sites become 
healthy at rate 1, and they infect simultaneously all their neighbors at rate A.) The 
behavior of the process is not expected then to be qualitatively different from that of the 
much studied (linear) contact process (see Chapter VI of [Ligl] and Part I of [Lig2]). For 
this reason we will focus in this paper on the cases ^ > 2, which are more challenging and 
do present a different qualitative behavior. As we will explain in the next subsection, 
this difference in behavior is indicated by the associated mean-field model. In the end of 
this introduction we will review some results about threshold 9 >2 and related models, 
from [Toom], [DC], [BG], [Chenl], [Chen2]. 

For coupling purposes, it is convenient to construct the process using a system of 
Poisson marks. For this purpose, associate to each site in V two independent Poisson 
processes: one with rate 1, and one with rate A. Mark the arrival times of the former 
with symbols D (for "down") and those of the latter with symbols U (for "up"). Make 
these Poisson processes independent from site to site. Use the marks now in the obvious 
way, to define the process: A spin 1 at site v fiips to when it encounters a D mark there; 
a spin at site v fiips to 1 when it encounters an U mark there and at least 9 neighbors 
of V have spin 1 at that time. The probability space on which these Poisson processes are 
defined will be large enough to accommodate the process started from arbitrary initial 
configurations. 

We will denote by {tJq q x.^t>o the process started from a random distribution picked 
according to law ^ at time 0. When jj. is product measure with density p we will use the 
notation (ry^ ^ When there is no risk of confusion, G, 9 and A may be omitted from 
the notation. 

The point mass on the configuration with all sites in state i e {0, 1} will be denoted 
5G,i- The distribution 5g,q is trivially invariant for the threshold 9 contact process, when 
^ > 0. By attractivity, 77^ q ^ ^G,e,\i as t — > 00, where =^ denotes convergence in 
distribution, and i>G,e,\ is called the upper invariant measure. 

We say that the process started from the distribution dies out when 77^ ^ ^.^ ^ ^G.o? 
as t — >^ 00. When this happens for every /i, we simply say that the process dies out. 
Attractivity implies that the process dies out precisely when i'G,e,\ — Sg,o- When the 
process does not die out, we will say that it survives. 



THRESHOLD CONTACT PROCESSES 



3 



For V eV, set 

We will use for ^ the same conventions on notation as for tJq g Also, when 

G is such that Pq q x-t^'^) does not depend on v, we will omit v from this notation. The 
critical point for the threshold 9 contact process on G is defined by 

Ac(G, 9) = sup{A : i'G,e,x = Sg,o} = sup{A : for each v eV, PG,6>,A;t(^) ^ as t ^ oo}. 

Clearly, the convergence of Pq q x.t{v) to cannot be faster then exponential. Explicitly: 

Pg e A;t(^) ^ P(there is no D mark at v from time to t) > e~*, (1.1) 
for t > 0. It is natural to define 

Aexp(G, 9) — sup{A : for each v eV ^ PG,e,X;ti'^) ~^ ^ exponentially fast as t — > oo}. 
Obviously 

Aexp(G, ^) < Xc{G,9), 

and it is interesting to decide when equality holds. 

Even when the process survives, it may happen that for small p > 0, 1]^ ^ ^.^ =^ Sg,o, 
as t — s> oo. By attractivity, if this happens for some value of p, it will also happen for 
smaller values of p. This leads to the following definitions: 

PciG, 9, A) = sup{p e [0, 1] : VG,e,\;t =^ ^g,o, as t ^ oo} 

= sup{p e [0, 1] : for each v e y, ^ j^.^(v) — > as t — > oo}, 

Pexp(G, 9, A) = sup{p e [0, 1] : for each v e V, pQ g x-ti"^) exponentially fast as t — > oo}. 
As above, obviously 

Pe.p{G,9,X) < Pc{G,9,X), 

and again it is interesting to decide when equality holds. 

When (rjQ q x-t) survives, it is natural to ask if it converges in distribution as t ^ oo. 
This is in general a difficult question, but the following provides a partial answer and, 
by-product, an estimate on Pc{G, 9, A). 

Proposition 1. For any G, 9 and \, for p > A/(A + 1), 

VG,e,X;t ^ ^G,e,x, as t ^ oo. (1.2) 
Therefore, if the process survives for a certain value of X, then 

PciG,9,X)<j^. (1.3) 
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Proof. Let /3g,p be the product measure with density p. For any ^ > 0, the threshold 
6 contact process is stochastically dominated by the threshold contact process. This 
latter process is simply an independent flip process whose unique invariant distribution 
is pG,x/iX+i)- Hence, for p>X/{X + l), 

T^G,e,x < /3g,a/(a+i) < I3g,p < 5g,i-, stochastically. 

But as t ^ oo, we know that VG^e x-t ^ ^G,e,x ^ind Vg e x t ^ ^G,e,x- Therefore (1.2) 
also holds. ■ 

Remark on case ^ = 1 .■ In that case, duality can easily be used to show that in Proposition 
1, (1.2) holds for every p > and hence (1.3) can be replaced by Pc{G, 9, A) = 0. 

1.2 The mean-field model. When G is regular of degree k (i.e., each site has k 
neighbors), its is natural to compare the evolution of the threshold contact process on 
G with a corresponding "mean-field" evolution. By this we mean the evolution of a 
deterministic density ip^^'x-t)t>o, which is governed by 



d 



dt 

with p]f^'x-o = Pj ciiid where 



MF,p MF,p I -I /-I MF.p \ -D- f MF,p n\ /i a\ 

Pn,e X;t = -P^,e X;t + ^ (^^ " P^/A;* j ^M.A;*' ' (l'^) 



is the probability that a binomial random variable with k attempts and probability of 
success X is larger than or equal to d. 

To explain the origin of this mean-field evolution and its relationship with the thresh- 
old contact process on G, we observe that from the definition of that process, for each 

j/G,e,xM = -p'gax-A^) +A P (v'gax-A^) = 0, # e A/^,. : v'g^x-A^) = l} > ^) • 

One then obtains (1.4) if one pretends that the k, + 1 random variables Vg e x-ti''^) ^ 
TjQ Q x-tiu), u E N'g,vi are i.i.d., with common density p^^'x-t- 

Trivially, pfJ'J-t = 0, for aU t > 0. When p e (0, 1], it clearly follows from (1.4) that 
p^^'x-t ^ P^~^ > 0' t>0. It will then be convenient to rewrite (1.4) as 



dt 
where 



'^iog(pSt) = ^('^'^'^;pSt)' (1-5) 



H{k, e, A; a;) = -1 + A^^^Bin(K, x, 0), x e (0, 1] 
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Note that when 9 > 2, then < Bin(/t, x, 6) < {k'^/2)x'^, and hence 

lim H{k, 9, X; x) = -1. (1.6) 

a;\0 

Also important are the elementary facts that H{k, 9, A; x) is continuous in x and in 
A, H{k,9,X;1) — —1, H{k,9,0;x) = —1, and, provided k > 9, for each x e (0,1), 
H{k,, 9, A; x) is strictly increasing in A, with lim>^^oo H{k, 9, A; x) — oo. 
These facts motivate the definition of the critical point 

A^^(k, 6*) = sup{A > : sup H{k, 9, A; x) < 0}, 

a:e(0,l] 

and the critical density 

P^^{k, 9, A) = inf{x G (0, 1] : H{k, 9, A; x) > 0}. 
Note that the facts above imply that, when k, > 9 >2, 

< X^^{k,9) < oo, 

and 

< p^^iK, 9,X)<1 for A > Xf^{K, 9). 
Moreover p^^{k, 9, A) is strictly decreasing in A > X^^{k, 9), with 

]imp^^{K,9,X) = 0. (1.7) 

A^oo 

Define 

V^^{k,9) = {{X,p)e[0,oo]x[0,l]:X<X^^{K,9)orp<pf^{K,9,X)}. 

Proposition 2. For every k>9>2, the following dichotomy holds. 
In case {X,p) e {V^^{k,9))^, 

^ Pc'''i^,0,X) > 0, forallt>0. (1.8) 

In case {X,p) e T>^^{kj 9), then for some C e (0, oo), 

pe-*< p^^'l^ < Ce-\ forallt>0. (1.9) 

In particular, a discontinuous transition happens at X^^{k,9). 

Proof. The only statement that requires explanation is the upper bound in (1.9). To 
prove it, note first that (1.5) implies that when (A,p) e T>^^{k,, 9) we have limt_>oo P^j q 'x-t - 
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0. Using (1.5) a second time now, this time in combination with (1.6), shows that for 
any e > 0, there is C e (0, oo) such that 

^ C"e-(i-^)*, foralH>0. 

Using this estimate in combination with (1.4) now yields, since 9 > 2, 

d MF,p ^ MF,p , X '^^ p-2(l-e)t 

Multiplying by e*, we obtain 

jMpZ%:) i A^(C')^e-<'-<)'. 

Supposing e < 1/2, integration in t from to s yields 

rs 2 /-oo 2 

e pK,e'x;s - P < A Y (C^ ) e ^ ' dt < X — {C) e ^ > dt = C < oo. 

(1.9) follows, with C = C" +p. ■ 

It is easy to see that A^^(re, 6*) — > 0, as — > oo. The rate at which this convergence 
occurs is also easily identified from standard facts about convergence of binomial distri- 
butions to Poisson distributions. For this purpose extend the definition of Bin(K, 6) 
to be 1 when a; > 1, and observe that straightforward computations then yield 



. ^ _ 7 

Bin(K, 7/k, 6*) — > Poisson(7, 6*) := cis — >• oo, (1-10) 



.7' 
i>e 

uniformly in 7 > 0. The corresponding limit for the function H{k, 9, A; x) is 

H(K,9,(j)/K;^/K) — > — 1 + (f) — '■ — — ^, as K — > 00, 

7 

uniformly in 7 > 0, for each (p > 0. From this it is easy to derive 

lim KXf^{K,9) = $e inf ^ . ^ , (1.11) 
«^oo c V w Poisson(7, 6*) 

The constants ^0 can easily be shown to satisfy 

1 < $2 < $3 < ... and lim ^0 = 00. 

e^oo 

Remark on case 9 = 1: In this case, in contrast to (1.6), we have lima;\^o H{k, 9, A; x) = 
— 1 + Xk,. It is an instructive exercise to use this fact and the bound H{k, 9, X; x) < 
—1 + Xk, for a; > 0, to analyze the behavior of the mean-field model in this case. In 
contrast to Proposition 2, one finds a continuous transition at A^^(k, 1) = 1/k, with 
P^^{k, 1, A) = for aU A > A^^(k, 1). The analogue of (1.11) is also true, with $1 = 1. 

1.3 Results for the process on homogeneous trees. We will study threshold contact 

processes on the homogeneous tree T;,, of degree 6+1. 

Our first result provides conditions for survival of the process on T;, based on the 
survival of the mean-field model with k = b (note: not k = b + 1) and A/ (A + 1) in place 
of A. 
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Theorem 1. For arbitrary b>2 and 9 >2, ifX^^{b, 9) < 1, then 

and for arbitrary A > X^^{b, 9) /{I- Xf^{b, 9)), 

Pc{T,,9,\) < p^^{b,9,X/{^+l)) < 1. (1.13) 

Since we know that lim^^oo A^^(6, 9) = 0, Theorem 1 contains meaningful statements 
when b is large. Moreover, combining this theorem with (1.11), we learn that 

limsup bXc{Tb,9) < ^q. (1.14) 

6— »oo 

This result will be sharpened in Theorem 3 below. In our approach to the proof of 
Theorem 3, we will prove first the somewhat technical Theorem 2 below. Note that 
thanks to (1-14), this theorem covers the behavior near Ac(Tb, 9), when b is large. This 
theorem should be compared to Proposition 2. 

Theorem 2. For arbitrary 9 > 2 and A G (0, oo), there are bo, 6 e (0, oo) such that if 
b > bo and A < A/b, then for every p G [0, 1] the following dichotomy holds. Either 

1 '■^ 



liminf — / pZ a^-^dt > -, (1-15) 
or, for some C G (0, oo), 

PnAX;t ^ ^e"'' forallt>0. (1.16) 

Moreover, for b > bo the set 

VA{Tb,9) = {(A,p) G [0,A/b] X [0,1] : alternative (1.16) holds} 

is an open subset of [0,A/6] x [0,1] in the relative topology induced by the Euclidean 
topology ofM?. 

Note that alternative (1.15) implies that the process survives when it starts with 
density p. For p= 1, (1.15) is equivalent to each one of the statements 

Pn,e,X;t ^ ^' foralH>0, 

and 

The following is immediate from Theorem 2, (1.14) and (1.3). 
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Corollary 1. For each 9 >2, the following statements hold when b is large. 

< Aexp(T6,^) = Ac(T6,^) < oo, 

and the process survives at this critical point. Moreover for A > Ac(Tb,6') close to this 
critical point, 

< Pexp(T6,^,A) = Pc{Tb,e,X) < 1, 
and the process started from this critical density Pc(Tfe, 9, A) survives. 
Note that in particular, under the conditions in Corollary 1, 

Pc{Tb,9,Xc{Tb,9)) < 1. 

Theorem 2 and Corollary 1 show qualitative similarities between the behavior of the 
threshold contact process on T5 and the corresponding mean-field model, when b is large 
and A is small. The next theorem shows a related quantitative similarity. 

Theorem 3. For arbitrary 9 >2, 

lim 6Aexp(T6,^) = lim6Ac(Tb,^) = ^0. (1.17) 

6— »oo 6— >oo 

Remark on case 9 = 1: In that case, duality can easily be used to show that (1.17) also 
holds, with $1 = 1. 

In contrast to the results above, for large values of A, the process on and the 
mean-field model behave differently, as the comparison between the following theorem 
and (1.7) shows. Set 

Pexp(Tb, ^, 00) = lim Pexp 

{Tb,9,X), p,{Tb,9,oc) = lim pe(Tb,^,A). 

A— >oo A— >oo 

Theorem 4. For arbitrary b >2 and 9 >2, 

Pc{Ti„9,(X)) > Pe^p{Ti„9,(X)) > 0. 

We do not know if the critical densities pd^b^ 9, 00) and Pexp(T6, 6*, 00) are identical 
to each other, but the next theorem shows that at least they display similar behavior as 
b — > 00. 

Theorem 5. For arbitrary 9 >2, 

< liminf 6^/^^-^^Pexp(T6,^,oo) < limsup 6^/^^"^) pc(T6, ^, 00) < 00. (1.18) 

6-»oo 6^00 

1.4 Results for the process on oriented homogeneous trees. As a tool, for 
comparison purposes and for its own sake, we will also study the threshold contact 

— * 

process on an oriented graph T5, obtained from T5 in a fashion described next. 
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First we introduce some notation. We embed a copy of Z in T5 and use the notation 
£ = £0 to denote the set of vertices of T5 covered by this embedding. Abusing notation, 
we will denote the elements of L by the names of the elements of 1 that they represent 
in this embedding. We will also refer to site as the root of T;,. Define, inductively in 
n'>\,Ln^^ the set of vertices which are neighbors to some vertex in Ln-\ and are not 
in VSlZlU. 

With this notation, include iy^u) in the set of oriented edges of if v,^ e £ and 
tt = u + 1, or if V e jCn-i cind u G Ln-, for some n > 1. 

Theorems 1 to 5 have analogues for the threshold Q contact process on T;,. Theorem 
1 for Tfe is actually a corollary to the same statement for T5. Theorem 2 and Corollary 
1 admit much stronger versions for those are stated as Theorem 6 and Corollary 2 
below. The analogues of Theorems 3, 4 and 5 for Tb are also true and can either be 
obtained from the corresponding results for T?,, or can more easily be proved directly. 

Theorem 6. For arbitrary 9 > 2 and h > 2, the following dichotomy holds for every 
A > and p e [0, 1]. Either 

Pf„e,x;t ^ ^'c''^(^'^'^)' foraUt>0, 
or, for some C G (0, 00), 

< Ce-\ forallt>0 

Moreover, for every 9 > 2 and b > 2, the set 

V{fb,9) = {{X,p) e [0,00) X [0,1] : alternative (1.20) holds} 

is an open subset of [0, 00) x [0, 1] in the relative topology induced by the Euclidean topology 
ofR\ 

Corollary 2. For each 9 > 2 and b > 2 for which Aexp(Tb,6') < 00, the following 
statements hold. 

Xf^ib,9) < AexpA,^) = Xcifb,9), 
and the process survives at this critical point. Moreover for A > Ac(Tb, 9), 

pf^{b,9,X) < Pexp(T6,^,A) = Pc{fb,9,X) < 1, 

and the process started from this critical density Pc(Tb, 6*, A) survives. 

We do not know if the first inequality in each display in Corollary 2 holds for the 
threshold contact process on T;,. But it is known that the second of these does not hold 
for the process on Z'^, d > 3, with 9 = 2, as reviewed in the next subsection. 

1.5 Related previous results. Threshold contact processes with 9 = 2 and closely 
related models have been studied in [Toom], [DG], [BG], [Chenl], [Chen2], sometimes 
under the name "sexual contact process", and mostly on Z^. In [Toom] discrete time 
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versions were studied, and contour arguments were used to show survival. In [DG] 
these contour methods were adapted to continuous time; their main result can be stated 
as follows using our terminology. Let 71? be the oriented graph obtained from 71? by 
setting A/"g2 ^ = {t; + (1, 0), + (0, 1)}. It is proved in [DG] that Ac(Z2, 2) < oo. Note 
that, by stochastic domination, this implies that Ac(Z'^,2) < cxd, for (i > 2. In [BG] a 
renormalization procedure was introduced, which can replace the contour methods in 
proving survival. In [Chenl] and [Chen2] continuous time models which can be seen as 
modified threshold d = 2 contact processes on were studied. In one of these modified 
models the flip rates at v e Z*^ at time t > are given by 

• 1 flips to at rate 1. 

• flips to 1 at rate A in case there are at least 2 sites of Mz^^v that are separated 
from each other by Euclidean distance and are in state 1 at time t. 

The most important result from [Chenl], [Chen2] in connection to the current paper 
is the fact that for this modified model in d > 3, when A is large, survival occurs 
starting from any positive density p. Since the threshold 2 contact process dominates 
that modified model, we learn that when (i > 3 and A is large, pc(Z'^,2, A) = 0. This 
means that the qualitative behavior of the 9 = 2 threshold contact process on Z*^, d > 3, 
deviates from that of the corresponding mean- field models (for which ^^^^(2^, 2, A) > 0, 
for all A > 0), but this deviation is in the "opposite direction" of the deviation observed 
in the corresponding models on homogeneous trees (for which, contrary to the mean-field 
model, limA^ooPc(T2d-i,2, A) > 0). 

The only results that we are aware of for threshold 9 > 2 contact processes on trees 
are the following ones, from [DG]. There the authors consider the model with 9 = 2 
on T2. They state that the contour methods used in that paper can be used to prove 
that this process survives when A is large. They then prove, using this result, that the 
transition at Ac(T2,2) is discontinuous (our Corollary 2 extends this result). 

1.6 Organization of the paper. In Section 2 we prove the results about the threshold 
contact process on Tb (obtaining Theorem 1 as a corollary). In Section 3 wc prove 
the results about the threshold contact process on T;, when b is large and A is small, 
namely. Theorems 2 and 3. In Section 4 we prove the results about the threshold contact 
process on when A is large, namely. Theorems 4 and 5; for this purpose the bootstrap 
percolation model will be introduced as a tool. 

— * 

2. Comparison between the model on T5 and the mean-field model. In this 
section we will prove the analogue of Theorem 1 for T;, and Theorem 6. Both theorems 
result from a fairly direct comparison with the mean-field model. This comparison is 
based on writing down, for an arbitrary site v, the differential equation 

. = -P- . +APfe , (v) = 0, ^lueXf :?7? , (u) = l}>9), 

(2.1)) 

and noticing that ¥5 has the special property that the random variables rK. (u), 
u e J\ff ^, are independent and have the same distribution as r/S (w). If the random 
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variable ^ ^ ^(^) were also independent of those, (2.1) would reduce to the mean-field 

equation (1.4), but this independence does not hold. In each one of the two proofs below 
we deal with this lack of independence in a different way. 

Proof of Theorem 1 and its analogue for T^. We will prove that for arbitrary h>2 and 
e>2, if P(fe, 6) < 1, and A > Xf^{b, d)/{l - \f^{h, 9)), then 

Pc{%,0,X) < p^^{b,e,x/{^+i)) < 1. (2.2) 

This suffices, since it obviously implies 

AcA,«) < T^5^ < (2.3) 

and (1.12) and (1.13) follow respectively from (2.3) and (2.2), since the threshold 6 
contact process on T;, stochastically dominates the threshold 9 contact process on 

Only the first inequality in (2.2) needs to be proved. For this purpose consider the 
arbitrary site v that appears in (2.1) and define the following event E: rjZ ^ ^ ^(v) = 0, 

and either there is no U nor D mark at v between times and t, or else, the last such 
mark is a D mark. By using the time reversibility of Poisson processes, a standard 
computation gives 

nE)>^{l-p). (2.4) 

The event E is clearly independent of the random variables Vfg^ t^^-^' ^ ^ "^fi v' ^i^*^^ 
also E C {ryS (v) = 0}, it follows from (2.1), (2.4) and the observation after (2.1) 
that 

It is convenient to rewrite (2.5) as 

^A.e.., ^ i(M,A,P;4^,,,,,,), (2.6) 

where 

L{h,e,\p]x) = xH{b,e,\/{\+l)]x) + Bin(6,x,6') (x -p), xe[0, 1]. 

When Xf^\h, ^) < 1 and A > \f^{h, e)/{l - \f^{h, 9)), then A/(A + 1) > A^^(6, 9). So 
p^^{b, 0, A/(A + 1)) < 1, and for p > p^^{b, 9, A/(A + 1)) arbitrarily close to 
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pf^{b, 9, A/(A + 1)) we have H{b, 6, A/(A + 1); p) > 0. It follows then that also 
L{b,d, X,p;p) > 0. We claim that 

Indeed, set 

If (2.7) were false, we would have tp < oo. Then by the continuity of in t, 

we would have „ . = P, and d/dtp% < 0, at t = tp. But (2.6) implies 

d/dtp^ > L{b,6, X,p;p) > 0, at t = tp. This contradiction proves (2.7), which 
implies 

Pe(T6,^, X)<p. 

Since p can be taken arbitrarily close to Pc^^{b, 9, A/ (A+1)), the proof of (2.2) is complete. 
■ 

Proof of Theorem 6. Applying Harris' inequality to (2.1), we obtain 

jA,.., ^ ^ bl,../^^ = O) P (# {« . K : ,l^^ Ju) = l}>0). 

From the observation after (2.1), now 

U. . < + A (l - 4^^,^^ J Bin ((,, 4^^,^^^, 0) , (2.8) 



which for p > is equivalent to 

i'^Pi.J S «(M,A;4^,,,J. (2.9) 

If p = 0, then (1.20) holds. Suppose that p > and (1.19) fails. Then (2.9) implies 
that p^ ^ ^ ^ — > as t — > oo. The proof that (1.20) holds then can be completed as the 

proof of Proposition 2. This shows that for each A and p either (1.19) or (1.20) holds. 

The statement about the set V{Ti,,9) follows now from the fact that the negation of 
(1.19) is a "finite-time condition": 

If (2.10) holds for some (A,p), then, by continuity, it also holds close to this point, with 
the same t. ■ 
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3. The regime of small A. In this section we will prove Theorems 2 and 3. We will 
abbreviate the notation, omitting T5, 6 and A for instance in: 

We will compare the threshold contact process on T5 with the similar process in which 
the spin of one of the neighbors of the root is frozen in the state 1. Recall the definition of 
£ from Subsection 1.4, and the corresponding terminology and notation. In our modified 
process the fiip rates are as in the threshold 9 contact process on T^, except for the site 
— 1, where the spin is kept frozen in the state 1. We will start this modified process with 
each site other than site —1 taking independently the value +1 with probability p, and 
being in state 0, otherwise. The notation r/^'^ will denote this process. Define also 

a'f = nvV^{l) = 1), 

and abbreviate = cr°'^. 

For comparison, we will also consider the trivial threshold contact process on T;,, 
i.e., the process in which the spin of each site flips independently of anything else, with 
flipping to 1 at rate A, and 1 flipping to at rate 1. Let 7rf be the probability that in 
this process a given site is in state 1 at time t, when at time this probability is set to 
ttq = p. It is elementary that for every p e [0, 1], 7rf < tt^ \ A/ (A + 1) as t — > 00. In 
particular, there is i{X) such that 

Trf < A, for t > i{X). (3.1) 

By attractivity, for any p e [0, 1], t > and < Zi < I2, 

Lemma 1. For arbitrary b > 2, > 2, p e [0,1] and t >0, 

^ai'" < -ai'P + Xbai+^'P, l>0 
dt * ~ * t ■> - 

And 

j^a'f < -a'f + X{,,lbal + Bin{b,ale)}, 1>1. (3.4) 

Proof. We will use the following Terminology. For each site v of T^, the b sites in ACs „ 
will be called forward neighbors of v and the single site in Afjf,,v\-^j^^ ^ will be called the 
backward neighbor of v. 

From the definition of (r/^'^) and cr|'^, 

^^al'^ = -a'f + A P {r^V^{l) = 0, # e M : vViu) = l] > 9) 

< + A P ( # {« e A/i : r^l^^iu) = 1} > ^) . (3.5) 



(3.2) 



(3.3) 
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Inequality (3.3) follows from (3.5) and the observation that, since 6 > 2, for the 
site I to have at least 6 occupied neighbors, it must have at least one occupied forward 
neighbor. 

To derive (3.4) from (3.5), we compare the process 1]^^ with a further modified process 
in which the spins at the sites —1 and / are both frozen in the state 1, while the spins at 
other sites evolve as in the threshold 9 contact process. In this modified process, the spins 
at the neighbors of I evolve independently of each other. Note that by attractivity, the 
distribution of {r]p^{u))ueAfi is therefore stochastically dominated by a product measure 
in which the forward neighbors of I have probability af of being occupied, while the 
inside neighbor of / has probability Trf of being occupied. The probability in the r.h.s. of 
(3.5) is now estimated from above by the probability that either the backward neighbor 
of I and at least one of its forward neighbors are both occupied (recall ^ > 2), or else 
that at least 9 of its forward neighbors are occupied. ■ 

Lemma 2. For arbitrary 9 > 2 and A e (0, oo), there are b*,S*,t* e (0, oo), such that 
if b > b* and A < A/b, then for every p e [0, 1] the following dichotomy holds. Either 

ol > —, forallt>t*, (3.6) 

or, for some C e (0,oo), 

erf < Ce-°-^*, forallt>0. (3.7) 



Remark: The exponential rate 0.6 in (3.7), could be replaced with any rate smaller 
than 1, with minor modifications in the proof and a larger value for b* . In our proof of 
Theorem 2, all that we will need about this rate is that it is larger than 1/2. 

Proof. For later convenience, we take b* = 9A^. We will use Lemma 1, and for this 
purpose we need to estimate yrf . Under the assumptions in the lemma that we are 
proving, we have A < A/b < A/b* — 1/(9A^) =: A. Let (7rf)t>o be defined in the same 
way as (7rf)t>o, but with A replacing A. Clearly Trf < tt^ . So, by (3.1), there is t* = i{X) 
which depends on A, but not on 6 or A (once they satisfy the conditions in the lemma), 
such that Trf < A = 1/(9^^), for every p e [0, 1] and t > t*. 

We use now the two inequalities in Lemma 1. The first one with I = and the second 
one with / = 1. We suppose that t > t*, so that we can use the estimate above on 7rf . 
Since also A < A/b, these inequalities read then 



d 
Jt 



Multiply the first of these inequalities by 1/{3V A) and the second one by y A, and add 
the resulting inequalities to obtain 

d 2 A^/"^ 

j^xt < --xt + -^Bm{b,af,9), t>t*, (3.8) 
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where 

By (3.2), al'^ < erf. Therefore we obtain the foUowing comparison between xt and erf: 



From (3.8), the fact that Bin(6, a^, 9) < h^{pXf (since 6* > 2) and the first inequahty 
in (3.9), 

j^xt < -\xt ^ \A''l''h{xtf =: G{xt), t>t*. 

Note that G{x) < for < a; < 4/(9^^/^6), and hm^^o G{x)/x = -2/3. One can use 
these facts, as in the proof of Proposition 2, to conclude that if there is some i > t* such 
that 

< 51^. (3.10) 

then there is some C < oo such that 

Xt < C e-^-^\ for all t > 0. (3.11) 
But from (3.9), the condition (3.10) is implied by 

at < 



* 3A2(1 + 3A)6' 

and (3.11) implies (3.7) with C = Sy/AC. This completes the proof of the lemma, with 
6* = 4:/ {^A^ {1 + 3 A)). ■ 

Lemma 3. For arbitrary 6 >2 and A e (0, oo), there is hi e (0, oo), such that ifb>bi 
and A < A/b, then for every p e [0, 1] and I > 0, 

P?>^;'"-(^J (3.12) 



Proof. To derive (3.12) we consider the discrepancies between the processes {rjf) and 
{rjf'^). We construct these two processes using the same structure of Poisson marks, and 
for every site v of T;, we set 

In this fashion, 5t takes the value 1 at the sites where the processes disagree at time t, 
and at the other sites. 

Observe that a D mark eliminates a discrepancy, while a U mark will possibly create 
a discrepancy at a site v only if at least one neighbor of v has a discrepancy at the time 
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of this mark. Therefore 5t < Ct ' > where {Q' ) is a process in which the spin at the site 
— 1 is frozen in the state 1 while all other spins are initially set to 0, and evolve with the 
following rules at each site v ^ —1: 

• 1 flips to at Poisson D marks at v. 

• flips to 1 at Poisson U marks at v if and only if at least one of the spins in A^^, 
is in state 1 at the time of that mark. 

Note that these flip rules are those of the threshold 1 contact process on Tf,, with infection 
parameter A. 

The threshold 1 contact process is stochastically dominated by the contact process, 
in which spins flip at rates: 

• 1 flips to at rate 1. 

• flips to 1 at rate A times the number of sites of My in state 1 at time t. 

We denote by {Ct ) the process on T;,, in which the spin at the site —1 is frozen in the 
state 1 while all other spins are initially set to 0, and then allowed to evolve according 
to the flip rates of this contact process. 

The chain of comparisons presented above implies that 



Let {if ) denote the contact process started from the conflguration in which only the 
site V is occupied (and no spin is frozen). Self-duality for the contact process implies 



The function u{-) has played an important role in the study of the contact process on 
Tfo. The proof of (3.12) will be complete once we argue that under our hypothesis, 



For this purpose we refer to results in Chapter 4 in Part I of [Lig2], where references 
to the original contributions can be found. The contact process on b > 2, has 
two critical points < Xi{b) < X2{b) < oo. Theorem 4.1 in Part I of [Lig2] tells us that 

Mb) > l/i^Vb)- Therefore we can flnd hi so that A < A/b < X2{b), when b > bi. Display 
(4.49) of Part I of [Lig2] tells us that u{l) < (/3(A))^ where /3(A) := limz^oo(w(/))i/'. 
Finally Theorem 4.65 in Part I of [Lig2] tells us that /3(A) < 1/Vb when A < X2{b). This 
completes the proof of (3.12). ■ 

Lemma 4. For arbitrary 9 > 2 and A G (0, oo), there are bo, S, S*,t* e (0, oo), such that 
if b > bo and X < A/b, then for every p e [0, 1] the following dichotomy holds. Either 
(3.6) and (1.15) both hold, or else (3.7) and (1.16) both hold. 




p\ = nm = 1) < nCt\i) = i) < nCt\i) = i). 



m:'\l) = 1) = P(ei'^(-l) = 1, for some s e [0,t]) 
< ¥{iP{-l) = 1, for some s > 0) 
= P(eW(/ + l) = l, for some s e [0,t]) 



: u{l + l). 
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Proof. Let b*, S* and t* be as in Lemma 2. We will take bo > b*, so that under the 
hypothesis of the lemma that we are proving we know from Lemma 2 that either (3.6) 
or (3.7) holds. 

Suppose first that (3.6) holds. Define 

pP = lim inf / pPdt, a^'P = lim inf ^ / a^'^dt, = a^^P. 

T^oo T Jo T-»oo T Jo 

From (3.3) and A < A/b, for each / > 0, 

A^l al+^'^dt > 1 / al^^dt + 



Hence, 

Aa^+^^P > a^^P. 
By induction in I and (3.6), we obtain now 

We use now Lemma 3, and for this suppose that b > max{6*, bi}. Then, from (3.12) and 
(3.13), we obtain 

- \VbJ - Alb \^bj - A-by S'\Vb 

Taking bo > max{6*, bi} large enough, we have A/Vb < 1/2, when b > bo. Hence there 
is I such that 

pP > 

~ 2Aib 

for all b > bo. We conclude that (1.15) holds then with 5 = 5*/{2A^). 
Suppose now that (3.7) holds. Prom the definition of (ryf ) and , 

j^pI = -pI +AP(r7r(0) = 0, ij^{ueHi:r^m = ^}>0) 

< -p\ +AP(#^eAA,:r7f(0) = l}>^). (3.14) 

We compare the process (r/f ) with the modified process in which the spin at the root is 
frozen in the state 1, while the spins at other sites evolve as in the threshold Q contact 
process. In this modified process, the spins at the neighbors of the root evolve inde- 
pendently of each other. Note that by attractivity, the distribution of {jft{u))u^MQ is 
therefore stochastically dominated by a product measure with density . Therefore, 
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for some C < oo, where in the last step we used (3.7). Multiplying both sides of this 
differential inequality by e* and integrating yields (1.16) (see the end of the proof of 
Proposition 2 for an identical estimate). ■ 

Proof of Theorem 2. Lemma 4 established the claimed dichotomy. It also implied that 
under the hypothesis of Theorem 2, (1.16) is equivalent to the negation of (3.6), i.e., 

a-f < for some t>t*. (3.15) 

The statement about the set ^^^(T;,, 6) follows then from the fact that if (3.15) holds for 
some (A,p), then it also holds close to this point (with the same t). (As in the proof of 
Theorem 6, this is a typical "finite-time condition" argument.) ■ 

Proof of Theorem 3. Since (1-14) has already been proved, we only have to prove that 

liminf6Aexp > ^e- (3.16) 

6— >oo 

For this purpose, let A <^e and \ = A/h. We will show that then 

P^oo < when 6 is large. (3.17) 

From Theorem 2 and the remarks after that theorem, we know that this implies that, 
when h is large, alternative (1.16) must hold and hence A < Aexp- Therefore 6Aexp > A, 
and since A can be taken arbitrarily close to (3.17) implies (3.16). 

From the proof of Proposition 1, in the introduction, we know that if p > A/(A + 1), 
then < p^, for every t > 0. Therefore (3.17) will follow once we show that 

inf < t:^, when h is large. (3.18) 
t>o ^ ~ 63/2 o ^ ^ 

To prove this claim, we use again (3.14), but this time we compare the process (ryf )t>o 
with the modified process in which the spin at the root is frozen in the state 0, while the 
spins at other sites evolve as in the threshold 9 contact process. We denote this modified 
process by (?7®'^)t>o- Let Ft be the event that the origin is vacant at time and that 
between time and time T there is no U mark at the origin. Then, for < t < T, 

P(#{«eA/i :r/f(«) = l}>^) < P(#{«eM:r/r(«) = l}>^, Ft) + P((Ft)'^) 

= ¥{#{ueNi: vf'^'iu) = l}>e, Ft) + P((Ft)^) 
< F{#{ueJ\fi: vf'^iu) = 1} > ^) + P + XT. 

In the modified process (?7®'^)t>o, the spins at the neighbors of the root evolve inde- 
pendently of each other. Note that by attractivity, the distribution of {'nf'^{'(J'))u£j^o is 
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therefore stochastically dominated by a product measure with density pf. Therefore, 
(3.14) yields 

j/t < -P't + A(Bin(6 + l,pf,^) + p + AT), < t < T. 

We will use this inequality with p = X = A/b and T = b^^^. We also change variables to 
Xt = {b-\- 1) pf*. The inequality above then implies 

d .^.(. xt A A 2A? ^ ,1/4 

^x. < -X. + ^Bm(6 + l,^,9j + y + ^, 0<t<6'/^ 

If (3.18) were false, then there would be arbitrarily large b for which 

a^t > 7T7^5 foralH>0. (3.19) 

b I 

We would then have, for the b for which (3.19) holds, 

- log(.,) < -1 + - Bin (6 + 1, ^, ») + ^ + ^, < t < ftV. (3.20) 

Thanks to (1.10), as 6 — > oo, the right hand side of this inequality, as a function of xt-, 
converges uniformly to 

Since A < $51, we have sup^>o-^(^) — ~^-> some C > 0. Therefore (3.20) yields, 
when b is large, 

dt ' - 2' - - 

This implies 

b 

which for large b contradicts (3.19). This contradiction proves (3.18), and completes the 
proof of (3.16). ■ 



4. The regime of large A. In this section we will prove Theorems 4 and 5. Our lower 
bounds on Pexp(T6,6', 00) will be obtained by comparison with a bootstrap percolation 
model that we describe next. 

The continuous time bootstrap percolation model on the graph or oriented graph 
G = (V, E) with threshold 9 and infection parameter A can be defined by taking the 
threshold contact process {'r]G,0,X;t)t>o on G, with same threshold $, and suppressing all 
the flips from 1 to 0. In other words, the bootstrap percolation process has flip rates at 
V & V at time t > given by 

• 1 flips to at rate 0. 

• flips to 1 at rate A in case there are at least 9 sites of Mg,v in state 1 at time 
t, and at rate otherwise. 
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We will denote by {(q q x-t)t>o ^^e resulting process, started from a random distribution 
picked according to law /j, at time 0. If one uses the same Poisson system of D and U 
marks to construct (rjQ q ^.^) and {C,q q y^), then, clearly 

V^,e,x;t ^ C^,e,x;v ^11 1 > 0. (4.1) 

It is also clear that (q g ^.^ is increasing in time, and therefore has a limit, q It 
is also clear that q a oo '^^^^ depend on A > 0, and that it can be obtained by the 
following iteration. Let Sq be the set of sites which at time are in state 1. Recursively 
define then 

yj{ve {Sn-iT : #{Mg,v n Sn-i) >e}, n> i. 

The sets Sn increase, and their limit is linSn = 'S'oo = e \-oo- (This iteration is often 
taken as the definition of bootstrap percolation in discrete time.) 

From (4.1), it follows that, for any A > 0, the sites that are vacant in (Sgo are vacant 
in the process {r]Q g ^^.J at all times. 

When G = Tb, the observation in the last paragraph can be used as follows. Consider 
the clusters of occupied sites in (^oo) i-c, the connected components of the subgraph of 
Th induced by the sites in Soo- It is easy to see that since d >2, the sites that belong to 
finite clusters of will eventually be in state in the process 77^ ^ Suppose that ^ 
is product measure with density p. If S^o contains almost surely only finite clusters, then 
for any A > the process {rjj^ q \-t)t>o '^^^^ Therefore p < PcC^b, 00). The next 
lemma provides a stronger conclusion under a stronger assumption. In its statement 
and its proof, we will use the following terminology and notation for T5. The distance 
between two sites is the length of the path that connects them. We will use for 
the ball of radius n and center at the origin. The outside neighbors of a site v are the 
neighbors of v that are farther apart from the origin than v (each site v ^ has b outside 
neighbors, and the origin has 6+1 outside neighbors). Denote by Rn the event that the 
site and some site separated from it by distance n are in the same cluster of Soo- 

Lemma 5. Suppose that for bootstrap percolation on T5, started from product measure 
with density p, P(i?n) decays exponentially with n. Then p < PexpCTfo, 9, 00). 

Proof. Let {rjj^g ^.^) be the threshold contact process started with the ball of radius n 
around the origin fully occupied and all other sites vacant. From the observations above, 

P (<,.,A;t(0) = 1, (Rnr) < P (<>,A;t(0) = l) • (4-2) 

The process {i^j^g ^.-j) is stochastically dominated by the process started from the 
same configuration, in which a spin never fiips and a spin 1 fiips to 0, at rate 1, iff 
all its outside neighbors are in state 0. For this process, let r^, f e Bn, be the random 
amount of time needed for the spin at v to flip to after the moment when it became 
allowed to flip. Clearly the Ty, v e B^, are i.i.d., with exponential distribution with 
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mean 1. A simple induction argument, starting from the sites at distance n from the 
root, and moving inwards, shows that the root will then flip to at the random time 



max y Ty, 



> t 



where is the set of (6 + 1) 6" ^ paths from to the sites that are at distance n from 
it. 

We obtain therefore, from (4.2), and the hypothesis of the lemma 

< Ci + P I max V 

\ v£-ir 

\vETf / 

where C\^Ci are positive finite constants and vr is an arbitrary element of n„. Taking 
n = \_et\ , for some e > small enough, a standard large deviation estimates for Poisson 
random variables (see, e.g., (A.l) in the Appendix of [KS]) shows that 



P J]t„ > t] < Cae- 



where C3 e (0, 00) and C4 is large enough that e~^-* < b. The last two displayed 
inequalities imply then that P ^77^^ q \ ti^) — decays exponentially with t, completing 
the proof. ■ 

Bootstrap percolation on homogeneous trees has been studied in [CLR] and [BPP]. 
Below we could build on some of their estimates. Nevertheless, for the reader's benefit, 
and at little extra cost, we will present a self-contained approach to our problem of 
estimating P(i?„), in order to use Lemma 5. 

To study the bootstrap percolation process on T;,, it is convenient to study also, as a 
tool, the bootstrap percolation processes on its subgraph T^, induced by the following 
subset of vertices. The set is the minimal set of vertices of T5 with the properties 
that e and if v e then Af^^ ^ C V^. We will also consider bootstrap percolation 

on the oriented graph T^, which has as set of vertices also , and defined then by 

First we observe that bootstrap percolation on and on are strongly related to 
each other in the following way. If we start them from a same set of occupied sites, Sq, 
then for n > 0, either both will have the root in Sn or neither one will have it. To see 
this, given a set of sites i? of T^, say that a site v E R is hidden from the root in R 
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if there is another site u e R which belongs to the path which connects the root to v. 
Observe that if in the iteration which defines Sn for either one of the two processes that 
we are considering we ehminate all the sites that are hidden from the root in 5^-1, we 
do not change the truth or falsehood of the statement that the root belongs to Sn- But 
with this modification, the sets Sn are the same for both processes. 

It is easy to write down a recursion for the probability p+ that the root belongs to 
Sn in the bootstrap percolation process on or T^, started from product distribution 
with density p. In the last paragraph we argued that is the same for both processes. 
Now, for the process on T^, the root will belong to 5"^ in case it belongs to Sq^ or in 
case it does not belong to So, but at least 9 of the sites in J\ff+ ^ are in Sn-i- This 

observation and some obvious facts about the geometry of yield: 

p+ ^ p + (l-p)Bin(5,p+_i,^), (4.3) 

with initial condition = p. 

The right-hand-side of (4.3) is an increasing continuous function of p^. Therefore 

Pn / Pto ■= mf{x > Q : X = p+ {1 - p)Bm{b,x,e)}. (4.4) 

The limit p'^ is the probability that the root belongs to 5*00 in this bootstrap percolation 
process on or T^. Since 6* > 2, p+ < 'mi{x > : x = p + {b/2)^x'^}. It is easy to 
use this observation to conclude that 

p+\0 asp\0. (4.5) 

The following concept will be used in the proof of Theorem 4. Consider bootstrap 
percolation onG = {V, E), and leiW <ZV and w e W. We will say that "tu is eventually 
VF-internally occupied" if w becomes eventually occupied in the bootstrap percolation 
process restricted to W. To make the definition precise, set Sq^ = SqH W, 

S^ = S^_,V}{veW\SZ, : #{A/G,.n5r-i}>^}, n>l. 

We now say that tv is eventually VF-internally occupied in case w E S^ := UnSY ■ Note 
that this event depends only on the initial configuration of occupied sites in W. 

Proof of Theorem 4- We will use Lemma 5, and for this purpose we need to estimate 
F{Rn). Consider bootstrap percolation on Tb and let {0 n} denote the event that the 
sites and n belong to the same cluster of Soq. Note that this is the same as the event 
that the sites 0, 1, n are all eventually occupied in this bootstrap percolation process. 

We will use the definition in the last paragraph before this proof in the case G = T5, 
W = C^. Note that the subgraph of induced by W (i.e., the subgraph of Tf, obtained 
by removing C from the set of vertices, along with the edges incident to these vertices) 
is an infinite collection of copies of . The roots of these copies of are neighbors to 
the sites in £, with each site in C being neighbor to 6 — 2 of these roots. For j e £, define 
Xj as the number of neighbors of the site j which are eventually VF-internally occupied. 



THRESHOLD CONTACT PROCESSES 



23 



When bootstrap percolation on Ti, is started from product measure with density p, it 
follows from the remarks above that each Xj has a binomial distribution corresponding 
to 6 — 2 attempts each with probability of success, where is given by (4.4). Clearly 
the Xj are also mutually independent. 

For each j E C define a grade as follows. If the site j is in state 1 at time 0, give 
this site grade A. If not, give this site the grade according to: also grade A if Xj > 9, 
grade B if X^- = ^ - 1, grade C\iXj = d- 2, grade YiiXj<d- 3. The probability of 
obtaining grades A, B or C are then, respectively: 



PA = P + (l-p)Bin(6-2,p+,0), 
PB = (l-p)Bin(6-2,p+,^-l), 

PC = (l-p)Bin(6-2,p+,^-2). (4.6) 

Observe that if {0 <-> n} occurs, then the following must happen: 

(i) No site in 0, 1, n can have grade F. 

(ii) If the sites j and k, with < j < k < n both have grade C, then there must exist 
a site i with j < i < k with grade A. 

For a given realization of the process, denote by ua, and nc, respectively, the number 
of sites in {0, n} which receive grades A, B, and C. Then (i) implies 

UA + riB + nc = n + 1, (4.7) 

while (ii) implies 

nA>nc- 1. (4.8) 

Since there are 3"'"'"^ ways to assign grades A, B and C to the sites in {0, n}, it follows 
that 

P(O^n) < 3"+^ max {pAT^iPBr^ {pcV 

nA,nB,nc 
nA+riB+nc—n+l 
nA>nc—i 

< 3"+^ max {pAT^{PBy{pcy 

nA,nB,nc 

< 3^^+^ max {pAy^iPBy 

nA,nB 
nA+nB>n/2 

< 3"+i (max{p^,j9B})"/'. (4.9) 
From the geometry of T5 and (4.9) we obtain 

F(i?n) < (6+1)6^-^0^71) < 3^^ [36(max{pA,M)'^^ 

L 

But (4.5) and (4.6) imply that for p > small enough, mai:x{pA,PB} < 1/(36)^, and 
hence 

P(i?„) — > 0, exponentially fast as n — > 00. 
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Lemma 5 now implies Pexp(Tb, 6*, oo) > p > 0. ■ 

We turn now to the proof of Theorem 5. The origin of the exponent 6 /{9 — 1) there 
is the following. Bin(6, 7/6", ^) is of order 1/6" for large h i& a — 6 / {9 — 1). The precise 
version of this statement that we need below is the following one, which can be checked 
by elementary computations. For arbitrary 7 > 0, 

m 

Z_^,e/(e-i)Bin(6,7/6^/(^-^),^) ^1, as 6 ^00, (4.10) 

uniformly in 7 G (0,7]. This has the following consequence for the mean-field model. 
For arbitrary 7 > 0, 

6*, A; 7/6^/(^-1)) ^ -1 + A as 6 ^00, 

uniformly in 7 G (0,7]. Therefore, for arbitrary A > 0, 

h'/^'-^^pf^{h,e,X) ^ {e\/Xf'^'-^\ as 6^ 00. (4.11) 
Proof of Theorem 5. From (1.13), in Theorem 1, 

Pc{Tb,e,oc) < ^^(6,^,1). 
Combined with (4.11), this implies 

limsup 6^/(^-1) Pc(T6,^, 00) < (^!)V(^-i), 
6^00 

which provides the upper bound in (1.18). 
The proof of the lower bound in (1.18), 

liminf6^/(^-iVe(T6,6',oo) > 0, (4.12) 

6— >oo 

builds on the proof of Theorem 4. We will use the same notation as in that proof. Define 
also 7 = p 6^/(^-1), 7^ = p+ 6^/(^-1), and 700 = pto 6^/^^"^^- The recursion (4.3) implies 

7n < 7 + b'/^'-'^Bm{b,-fr^-i/b'/^'-'\e) (4.13) 

Set 

7 = inf{a; > : a; = 7 + 2x^91}. 

(In this definition the factor 2 is arbitrary; any number larger than 1 could be used 
instead.) Since 9 >2, for small 7 > we have < 7 < 00, and 

(7)^ 

7 = 7 + 2^. (4.14) 
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Moreover, similarly to (4.5), 

7 \ as 7 \ 0. (4.15) 
Our next goal is to prove that when b is large 



From (4.10), there is b such that for < y < 7 and b > b, 



b'/(^-^)Bm{b,y/b'/('-'\e) < (4.17) 

Since 7 + 2x^ /6\ is increasing in x, when y < 7, we have 7 + 2y^ /9\ < 7. Note that 
7 < 7. Hence (4.13) and (4.17) imply, by induction on n, that 7„ < 7. Therefore 700 < 7 
and (4.16) follows. 

Combining (4.6), (4.16) and (4.15) we obtain, when b is large. 



6^/(^-1) 

IB 



PB < ^, 

1c 

PC < ^(g_2)/(g-i) : (4-18) 

where lim-y\^o7A = 111117x075 = hm-ys^o 7c = 0. Note that (4.18) implies the following 
technical estimate: 

7' 

max{{pAf,{pBf,PcPA}< ^, (4.19) 
where lim^x^o 7' = 0. This estimate is useful in combination with the following one: 

P(O^n) < 3"+i {ma^{{pA)MpBf,PcPA})^^^^^~'. (4.20) 

To prove this inequality, one can match pairs of sites in {0, n} in the following way. 
Recall that if the event {0 <-> n} happens, then the facts (i) and (ii) in the proof of 
Theorem 4 must happen. Prom fact (ii) wc know that each site which receives a grade 
C is followed eventually by a site with grade A, except possibly for the last site with 
grade C. Pair each site with a grade C with the first site with grade A after it, leaving 
possibly one unmatched site with grade C. Considering the sites with grade A which 
are unmatched to any site with grade C, we match the first of these sites to the the 
second one, the third to the fourth, etc, leaving at most one unmatched site with grade 

A. Pinally we match the first site with grade B to the second such site, the third site 
with grade B to the fourth such site, etc, leaving at most one unmatched site with grade 

B. Since the number of sites in {0, ...,n} is n + 1 and there are at most 3 unmatched 
sites, the number of matched pairs is at least ((n + 1) — 3)/2 = (n/2) — 1. The estimate 
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(4.20) now follows from the fact that the number of ways to assign grades A, B and C 
to the sites in {0, ...,n}, is 3"+^. 

Combining (4.19) with (4.20), we obtain, when b is large, 

HRn) < (6+l)6^-^P(0 ^ n) < (6+1)6^-^3"+^ ( = 27(6+1)6 {3^/y) 

By taking 7 > sufficiently small, we can make < 1. Then 

P(i?n) 0, exponentially fast as n — > 00. 
Lemma 5 now implies 6^/(^-1) Pexp(T&, 0, 00) > 6^/(^-1) p = 7 > 0. This proves (4.12). ■ 
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